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ABSTRACT 



The persistent Be/X-ray binary RX J0440.9+4431 flared in 2010 and 2011 and has been followed by various X-ray facilities (Swift, 
RXTE, XMM-Newton, and INTEGRAL). We studied the source timing and spectral properties as a function of its X-ray luminosity 
to investigate the transition from normal to flaring activity and the dynamical properties of the system. We have determined the orbital 
period from the long-term Swift/BAT light curve, but our determinations of the spin-period are not precise enough to constrain any 
orbital solution. The source spectrum can always be described by a bulk-motion Comptonization model of black body seed photons 
attenuated by a moderate photoelectric absorption. At the highest luminosity, we measured a curvature of the spectrum, which we 
attribute to a significant contribution of the radiation pressure in the accretion process. This allows us to estimate that the transition 
from a bulk- motion-dominated flow to a radiatively dominated one happens at a luminosity of ~ 2 x 10 36 ergs~'. The luminosity 
dependency of the size of the black body emission region is found to be tbb k L^' 3 02 . This suggests that either matter accreting 
onto the neutron star hosted in RX J0440.9+4431 penetrates through closed magnetic field lines at the border of the compact object 
magnetosphere or that the structure of the neutron star magnetic field is more complicated than a simple dipole close to the surface. 

Keywords. X-rays: binaries. X-rays: individuals RX J0440. 9+4431 . stars: neutron. 
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1. Introduction 

LS V +44 17/RX J0440.9+4431 was classified as one of the rare 
persistent neutron star Be/X-ray binaries based on the outcomes 
of the first ded i cated RXTE monitoring campaign performed by 
iReig & Rochel (l999). The source distance was later estimated 
at (-3.3+0.5) kpc bv IReig et alJ (l2005h . RX J0440.9+4431 is 
also reporte d in the seven-year and nine- year INTEGRAL all- 
sky survey (iKrivonos et al.l 1201 Obi |2012|) . From the former to 
the latter catalog, the average flux of the source increased from 
(0.95 ± 0.15)mCrab to (4.0 ± 0.1)mCratQ(17-60keV) as it un- 
derwent several episodes of enhanced X-ray activity in 2010 
and 2011 (see later in this section). The flux estimate given in 
the seven-year cata log is consistent with what has been mea- 
sured previously by Rei g & Rochel (1 19991) and can thus be con- 
sidered as the persistent level of the source. A similar flux 
(0.9 ± O.5)mCrat0 was also reported in the 54-month Palermo 
BAT Survey in the 15-150keV band (ICusumano et al.ll2010l) . 

The first observed outburst from RX J0440.9+4431, which 
we are awar e of, was detected on 2010 March 26, and lasted un - 
til April 15 dMorii et al.ll2010i iFinger & Camero-Arranzll2010l) . 
The Swift/BAT light curve (Fig. Q} showed that the source 
reached a peak luminosity of about 200mCrab (15-150keV) 
six days after the onset of the event. It remained in this state 
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1 In the energy band 17-60 keV these correspond to 1.35 and 5.7 x 
10~ 11 ergcirr 2 s -1 , respectively. 

2 In the energy band 15-150 keV this corresponds to 2.0 x 
10~" ergcnT 2 s -1 . 



for eight days, and then turned off to its persistent luminosity 
level in about ten days. The outburst was also observed by the 
RXTE All Sky Monitor (2-10 keV, ASM). In the ASM light 
curve, the source reached about lOcts/s (i.e., ~130mCrab). A 
detailed stu dy of the first outb urst from RX J0440.9+4431 was 
presented bv lUsui et al.ld2012l) . The authors show that the broad- 
band spectrum of the source could be described well by a com- 
bination of a black body component, a cut-off power law, and 
an iron line with a centroid energy of ~6.4keV. They also find 
a pronounced dip in the pulse profile interpreted as being due to 
obscuration of the X-ray emission from the neutron star by the 
accretion stream. 

A second outburst from RX J0440.9+4431 was detected by 
INTEGRAL on 2010 September 1 (Fig. [I). The measured aver- 
age and peak fluxes during the event we re ~30mCrab and (77 ± 
17)mCrab (17-60 keV), respectively (IKrivonos et all l2010al) . 
The Swift/BAT and RXTE/ ASM light curves showed that the 
outburst had a rise and decay time of about five days, with a 
compatible peak luminosity of ~80mCrab. The analysis of the 
combined Swift, INTEGRAL, and RXTE data led to identify- 
ing a possible cyclotron absorption feature at ~30keV, suggest- 
ing a surface magnetic field of B 3 x 10 12 G for the accreting 
neutron star hosted in this system. This magnetic field intensity 
would be compatible with the measure d properties of the sourc e 
broad-band noise in the X-ray domain (Tsygankov e t al.ll2012l) . 

The two intense X-ray emission episodes recorded from 
RXJ0440. 9+4431 showed the typical properties of the so-called 
"Type I" outbursts, usually displayed when the neutron star in a 
Be/X-ray binary approaches the system periastron and interacts 



Article number, page 1 of 1 1 1 1 



with the companion equatorial disk. The presence of a neutron 
star in RX J0440.9+4431 is secured by the detection of pulsa- 
tions in its persistent and outburst emiss ion. The measu red pulse 
period increased from »202s in 1999 dReig & Rochelll999l) to 
«206 s during the latest observations, thus evidencing an average 
spin down trend. The first estimate of the source orbital period 
was reported by Reig ( 20 1 1 ) at ~ 1 50 d using the orbit-to-spin pe- 
riod relation for Be/X-ray binaries. This period would be com- 
patibl e with the recurrence time of the source outbursts (150- 
160 dlMorii et alJl2010tlkrivonos et alJl2010aHTsvgankov et alJ 
l20l2h . 

A deep (20 ks) XMM-Newton observation of 
RXJ0440.9+4431 was carried out on 2011 March 18 to study 
the p roperties of its persistent emission (lLa Palombara et all 
l20l2h . During this observation, the source displayed an average 
flux of 6xl0~ n ergcm _2 s _I (0.3-12 keV) and the correspond- 
ing spectrum could be well described using a power law plus 
a black body model. The estimated radius for the thermal 
emission component was compatible with the size of a hot-spot 
on the surface of a strongly magnetized neutron star (~10 12 G). 
Pulsations at «205 s were detected throughout the observation. 
No significant spectral variability with the spin phase could be 
measured, although an indication of hardening was found at the 
pulse minimum. 

In this paper, we report on the last detected outburst from the 
source which occurred on February 201 1 (Fig. [T}. On this occa- 
sion, we triggered our pre-approved RXTE target of opportunity 
observations and performed a simultaneous monitoring cam- 
paign at the softer X-ray energies with Swift/XRT to study the 
evolution of the source spectrum with luminosity and improve 
the measurement of the source orbital period. We have also re- 
analyzed all available RXTE Swift/XRT, XMM-Newton, and 
INTE GRAL (iBradt et al.ll99l lGehrel s et al.l2004t l Jansen et al.1 
1200 it IWinkler et al.l 120031) observations performed in the di- 
rection of RXJ0440.9+4431 in 2010 and 2011 to compare the 
source emission properties across two orders of magnitude in lu- 
minosity. The data analysis technique is described in Sect. |2 
and all the results are reported in Sect. [3] Our discussions and 
conclusions are given in Sect.|4]and Sect. [5] respectively. 



2. Data analysis 

We have analyzed RXTE and Swift data using the heasoft soft- 
ware package (v. 6.12) while all s pectral fits hav e been per- 
formed using XSPEC version 12.6.0u (lArnaudll 19961) . 

For RXTE, we have followed the suggest ed analysis thread jp 
and we considered only PCA (2-60 keV, iJahoda et aLH l996) 
data from the PCU2 unit, which was always active. Good Time 
Intervals (GTI) were created by imposing that at least 30 min- 
utes passed after the satellite exited the SAA, and the elevation 
on the Earth limb was >10°. We additionally imposed an elec- 
tron ratio <0. 15 and modeled the background using the provided 
estimator v3.8. Events were extracted in the mode GoodXenon 
(time resolution 0.95 /ys and 255 energy channels). The instru- 
ment specific response matrices were generated taking into ac- 
count the corrected energy bounds and the gain corrections per- 
formed by the EDS. The average PCA spectrum was obtained 
from the standard2 mode without further rebinning and fit- 
ted above standard channel 4 (~3keV) up to the higher energy 
at which the source signal was detectable (i.e., about 20keV, 
for short observations or faint source states, higher otherwise). 
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Fig. 1. Light curve of RX J0440.9+4431 and times of the observa- 
tions reported in this paper. Black points with errors: Swift/BAT daily 
averaged light curve. Red vertical bars: times of the RXTE observa- 
tions. Blue horizontal bar: time span of the INTEGRAL data used in 
this paper. Green diamonds: times of the Swift observations used in 
our analysis. The first vertical yellow bar indicates the occurrence of 
the first maximum in the light curve; the others correspond to the times 
of the following maxima predicted by our orbital period determination. 
MJD 55 280, 55 430, and 55 580 correspond to 2010, March 25, 2010, 
August 22, and 201 1, January 19, respectively. 



Spectra were assumed to have a systematic uncertainty of 0.5%, 
as recommended by the instrument teanfl 

For Swi ft/XRT (0.5-1 1 keV) , we followed standard analysis 
procedures dBurrows et al.l 120051) and considered both Window 
Timing (WT) and Photon Counting (PC) data, but limited our 
spectral range to 0.8-10keV to avoid calibration uncertainties. 
XRT spectra were grouped to have at least 50 counts in each 
energy channel using grppha. We extracted separated spectra 
for data accumulated in each satellite orbit and used the proper 
exposure map to account for bad columns, hot pixels, vignetting, 
and PSF losses. 

We extracted the sou rce spectrum from the INTEGRAL 
IBIS/IS GRI (20-100 keV , lUbertini et al.1l2003l) and JEM-X (5- 
30 keV. iLund et al.ll2003b data of the 2010 outburst taken from 
2010, September 1 to September 10 using th e Off-line Science 
Softw are (OS A, v. 10) distributed by the ISDC (ICourvoisier et al.l 
2003). The JEM-X spectrum is extracted in the standard 16 bins 
defined in OSA, while the IBIS/ISGRI one in user defined 62 
bins. 

XMM-Newton data have been processed using XMM Sci- 
ence Analysis System (SAS) vers ion 12.0.1. We hav e extracted 
only the EPIC-PN (0.5-10keV. IStriider et all l200ll) spectrum 
from a region of 0.6 arcmin centered on the source. Background 
and p ile-up corrections were neglected (lLa Palombara et aD 
120121) . The spectrum was adaptively re-binned using the SAS 
tool specgroup in order to have at least 25 counts per energy 
bin and, at the same time, to prevent oversampling of the energy 
resolution by more than a factor of three. 



3. Results 

3.1. Orbital period 

To search for the source orbital period, we exploited the daily av- 
eraged light curve obtained from the hard X-ray transient mon- 
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itor page of Swift/BAT^ These data cover almost continuously 
the time interval MJD 53414-56282 (2005 February 13-2012 
December 31), owing to the large field of view of the instru- 
ment and the spacecraft pointing strategy. We first used the 
Lomb-Scargle periodogram method in the range 32-200 days. 
The highest significance period is found at ~150 days. We have 
also exploited the epoch folding technique on the same data by 
assuming 32 phase bins and assigning weights of 1 /<x 2 to the in- 
dividual bins (here <x is the reported uncertainty on the source 
count rate). The most prominent peak (x 2 = 942) is found also 
in this case at ~150 days with evident harmonics at multiple fre- 
quencies (see inset in Fig.|2|. We have refined our determination 
of the orbital period by limiting our search with the above meth- 
ods around the most significant peak and using 130 phase bins. 
The best value is found at (150.0 ± 0.2) days. The uncertainty at 
68% c.l. on this measurement is obtained by bootstrapping 1000 
light curves with count rates randomly distributed within the un- 
certainty of the measurements, and then computing the standard 
deviation of the period determinations. We have folded the light 
curve at this period and obtained the profile shown in Fig. [2] 
This is reminiscent of what is usually observed from Be/X-ray 
binaries, with most of the emission concentrated during the pe- 
riastron passage (phase 0.6 in the figure). We have fit the folded 
light curve with a constant plus a Gaussian to estimate the time 
of the light curve maximum. The Gaussian is centered at orbital 
phase 0.641 ± 0.001, which can be converted into the time of oc- 
currence of the first flare by taking into account the best orbital 
period: it corresponds to MJD 55 293.2 + 0.3, i.e., the time of the 
brightest flare observed so far from RX J0440.9 +443 1 . We have 
plot this time interval in Fig.[T| left section, as a vertical shaded 
region. We have then propagated the uncertainty on the orbital 
period to the following two flares observed in our campaign and 
plot them as shaded regions in the other two sections of the light 
curve. 

Interestingly, when limiting ourselves to the time range be- 
fore MJD 55 280 (2010 March 25), the orbital modulation was 
not detectable and no prominent peaks are clearly visible over 
the red noise. This indicates that the X-ray active phase of the 
source began only in April 2010 when it was detected in out- 
burst for the first time, since the availability of modern X-ray 
monitors. 

3.2. Timing analysis and spin period measurements 

To search for the ~ 200 s source spin period and its possible vari- 
ations, we performed a timing analysis of all RXTE data (event 
mode GoodXenon). Swift pointings were on average too short 
to perform a similar analysis. We first extracted the background 
corrected PCA light curves in the 3-35 keV energy range (time 
bin 2 s), and barycenter-c orrected them using the known optical 
position (Reig et al. 2005j). For each observation we removed 
the variability on time scales longer than the expected pulsed 
signaQ and determined the pulse period with the epoch folding 
technique. 

We obtained the best determined period for all the avail- 
able RXTE observations using an optimal number of phase bins 
(from 16 to 64 depending on the S/N), and estimated the un- 
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6 To achieve this, we computed the sliding average of the light curve 
on intervals of 400 s every 200 s and then interpolated these values at 
the central time of each time bin to subtract them from the light curve. 
The resulting time series was then used as input to the epoch folding 
algorithm. 
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Fig. 2. Folded Swift/BAT light curve (15-150 keV) of 

RXJ0440.9+4431 at the best determined orbital period, (150.0 ± 
0.2) days. The inset shows the periodogram obtained from the epoch 
folding technique applied to the Swift/BAT light curve. 



Table 1. Measurements of the spin period of RX J0440.9+4431 ob- 
tained from all the RXTE observations analyzed in the present work. 



Time 
MJD 


Ar a 

ks 


p 

1 spin 
S 


1 <x error 

s 


55292.35 


1.11 


205.10 


0.14 


55298.23 


1.24 


205.03 


0.13 


55301.23 


3.51 


204.84 


0.13 


55447.02 


3.60 


204.89 


0.17 


55447.95 


2.61 


204.24 


0.06 


55449.06 


3.32 


205.51 


0.07 


55450.24 


2.60 


204.60 


0.11 


55451.75 


3.67 


205.38 


1.19 


55452.99 


3.66 


205.26 


1.36 


55594.16 


3.67 


204.65 


0.29 


55595.65 


3.66 


203.97 


0.26 


55596.50 


2.79 


205.42 


0.19 


55597.55 


3.07 


204.57 


0.25 


55598.45 


2.95 


205.00 


0.18 


55599.50 


2.81 


205.99 


0.20 


55600.55 


3.35 


205.78 


0.25 


55601.39 


3.59 


205.86 


0.18 


55602.37 


3.56 


206.29 


0.18 


55603.49 


2.83 


206.21 


0.24 


55604.34 


0.97 


205.25 


0.86 


55605.45 


2.45 


206.08 


0.45 


55606.49 


1.48 


206.01 


0.47 



Notes. < a > AT = (r cnd - T st ,„) 12. 



certainties by running 100 simulations of the input light curves 
with rates estimated from Poisson distributions centered on the 
real measured values. The results are reported in Table [1] Some 
variability in the measured values is evident from the table, but 
the limited number and accuracy of our determinations are not 
sufficient to disentangle unambiguously the effects of the orbital 
motion and/or possible accretion torques. 
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3.3. Broad-band spectral analysis 

The availability of quasi-simultaneous INTEGRAL, RXTE, and 
Swift data permitted to investigate the broad-band spectral prop- 
erties of RXJ04 40.9+4431 (note that part of the date we re al- 
ready presented: Tsvgankov et al.ll2012t lUsui et al1l2012l) . The 
broad-band quasi-simultaneous data set is shown in Table[2]wifh 
a conventional ID of the corresponding multi-instrument obser- 
vations: A, S, 1, 2, 3, 4. For the INTEGRAL data reduction 
we took advantage of the most updated IBIS/ISGRI calibrations 
included in the OS A v.icfl 

To compare our fin dings with those reported previously by 
Tsvg ankov et al.l d2012l) . and test the effect of the improved IN- 
TEGRAL calibrations, we used the same data set as in the above 
paper (INTEGRAL data from MJD 55440-55449, RXTE/PCA 
data from OBSID 95418-01-03-00 and Swift/XRT pointing OB- 
SID 00031690003, i.e. observation S in Table EJ. We restricted 
the selection of pointings to a maximum off-axis angle of 7° 
to minimize the systematic uncertainties on the IBIS instru- 
mental response. The total effective exposure time of the final 
IBIS/ISGRI spectrum was 56 ks. This corresponds to an ef- 
fective exposure time for JEM-X2 of 64 ks (owing to the lower 
dead-time of the detector). 

A spectral model comprising a cut-off power law and 
a black body component could fit the data reasonably well 
(X 2 /d.o.f. = 1.09/558). Normalization constants (fixed to 1 for 
Swift/XRT) were also included in the fit to account for source 
variability and instrument inter-calibrations. Photoelectric ab- 
sorption is a ccounted for using t he m odel phabs with abun- 
dances from lAnders & Grevessel (119891) . Some residuals are 
present at ~ 30keV, as reported by Tsv gankov et al.l (120121) . 
However, we found that an additional Lorentzian absorption fea- 
ture with <x=6 keV at ~30keV, as proposed by these authors, did 
not significantly improve the fit (y 2 /d.o.f.= 1.05/556, 33% chance 
improvement probability). 

A better description of the spectrum could be obtained 
by using the using t he general Comptonization model bmc 
dTitarchuk et al.1 fl997h with an exponential cutoff at high en- 
ergy (bmc x exp(-£ , /£ , F O id),^ /2 /d.o.f.=0. 92/559). Because of the 
general shape of the convolutional kernel, this model can be ap- 
plied to Comptonization spectra of both optically thin and opti- 
cally thick regimes, and considering the cases in which thermal 
or the rmal plus bulk Comptonization dominate (Farinel li et alj 
120131) . The emerging spectrum of the bmc model can be ex- 
pressed as 

Nbmc A5 
1 + A 

where G(E, Eq) is the Green's function of the Componization en- 
ergy operator, and S (E) is the seed photon black body spectrum. 
The quantity A /(A +1) empirically provides the fraction of seed 
photons which are Compton scattered; it is related to the geomet- 
rical configuration of the system and to the seed photon spatial 
distribution in the Comptonization region. Actually, log 10 (A) is 
the model parameter returned by XSPEC. The Green's function 
G(E, Eq) is a broken powerlaw oc E a+i for E < Eq and oc E~ a 
for E > Eq. The value of the powerlaw energy index a is related 
to the efficiency of the Comptonization process. From the nor- 
malization of the model, it is possible to infer the effective radius 
of the region in which the seed black body radiation is generated 
by using the relation Nbmc—L39/D^ q (here L39 is the luminosity 
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Energy [keV] 

Fig. 3. Broad-band spectrum of the quasi-simultaneous IN- 

TEGRAL/ISGRI (green), INTEGRAL /JEM-X2 (blue), RXTE/PCA 
(red), and Swift/XRT (black) observations performed in September 
2010. The best fit is obtained with the absorbed bmc with exponential 
cut-off model. The residuals from this fit are also shown in the bottom 
panel. 



in units of 10 erg/s and D10 the distance in units of lOkpc): 

Here C=l for a spherical geometry and C=0.25 for a circular 
slab. Assuming for RX J0440.9+4431 a distance of 3.3 kpc, the 

relation above can be written as rBB.km = -^844 Nbmc / C(kT)^ eY . 

The best fit parameters obtained using the absorbed bmc with 
exponential cut-off model are reported in Table [3] We found 
that this description of the source continuum emission permit- 
ted to constrain the optical depth of the previously suggested 
Lorentzian shaped absorption feature to <0. 1 (see the residuals 
from the best fit in in Fig. 0. 

To test further the applicability of the spectral model in- 
troduced above to RXJ0440. 9+4431, we carried out a sec- 
ond fit to the quasi-simultaneous Swift/XRT and RXTE/PCA 
observations performed in April 2010 (observation A in Ta- 
ble |5J. At this epoch, the source was brighter (F2_30kev - 
3.8xl0~ 9 erg s _1 cm -2 ) than in observation S, PCA data could be 
used from 3 to 45 keV and an iron line was needed to fit the data 
(we modeled this a s a thin addition al Gaussian component cen- 
tered at ~6.4keV: IUsui et a"Dl2012l) . No evidence was found for 
the Lorentzian shaped absorpti on feature (width 6keV, centroid 
energy <42keV) suggested bv lTsvgankov etail f2M% . We ob- 
tained a 90% c.l. upper lim it of 0.03 on its optical depth. We also 
note that lUsui et al.l d2012l) did not find evidence of any absorp- 
tion feature in the combined RXTE PCA and HEXTE spectrum 
of the source. We conclud e that the absorption feature reported 
by iTsvgankov et ail was probably due to systematic bi- 

ases in the joint analysis of data from several instruments and/or 
a inadequate modeling of the continuum emission. 

3.4. Spectral variability 

To study possible changes in the spectral parameters of the 
source at different luminosity states, we fit here all the avail- 
able RXTE/PCA data, which provided an unbiased monitoring 



Article number, page 4 of 1 1 1 1 



C. Ferrigno et al.: Outbursts of RX J0440.9+4431 
Table 2. Log of observations of the 2010-201 1 outbursts used in this paper. 



ID a 


OBS ID b 


Exp. 


Start 


Stop 


log 10 (Flux) c 


Flux c 








ks 


MJD 


MJD 




xl0~ 9 


XlO 35 




Q541 8-01 -01 -00 


3 1 


55?9? 3306 


55?Q? ^778 


-8 941 ?+0 0007 


5 738 

J. / JO 


74 766 




95418-01-02-00 


3.3 


55298.2106 


55298.2554 


-8.4228±0.0008 


3.778 


49.222 


A 


hhai i AQf\f\f\ i 
UUUjloVUUUl 


A Q 
4.0 




CCQAA f\A OA 

JJJUU.U424 










95418-01-02-01 


3.2 


55301.2178 


55301.2611 


-8.7541+0.0013 


1.762 


22.954 




95418-01-03-01 


1.0 


55447.0197 


55447.0311 


-8.738±0.002 


1.827 


23.804 




i£JS 


56 


55441.2576 


55448.1785 


- 


- 


- 


s 


JEM-X 


64 


55441.2576 


55448.1785 


-8.78+0.01 


1.66 


21.6 




00031690003 


5.7 


55443.0294 


55443.2021 


— 


- 


- 




95418-01-03-00 


2.8 


55447.9302 


55447.9728 


-8.894+0.002 


1.275 


16.615 




95418-01-04-00 


2.5 


55449.0456 


55449.0763 


-9.021+0.003 


0.953 


12.417 




95418-01-04-01 


1.5 


55450.2285 


55450.2463 


-9.316+0.006 


0.483 


6.296 




yj4io-Ui-U4-U2 


i . i 


Jj4ji. /JOJ 


Jj4ji. / JO J 


-iU.j /+U.UO 


f\ f\A 


U.jj/ 




Q541 8-01 -04-03 


1 


5545? Q8^0 


5545? QQfi^ 


-10 47+0 07 


034 


0.444 




7UJUT u 1 W 1 V-/V7 


1 ^ 

1 .J 


14Q1 

JJJ7T. It-7 1 


SSSQA 1680 

JJJ7T. 1UOW 


Q Q^7+0 00^ 


1 1 


1 S 0S1 




a/co /i a 1 a 1 a 1 


1 c 

1.3 


crcnc /Z A A A 

jjjVj.o4ww 


CCCAC /CCT/I 
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96364-01-04-00 


3.2 


55598.4385 


55598.4848 


-10.38+0.02 


0.042 


0.547 




96364-01-05-00 


2.9 


55599.4872 


55599.5304 


-10.41+0.04 


0.039 


0.504 




96364-01-06-00 


2.8 


55600.5345 


55600.5763 


-9.362+0.004 


0.434 


5.659 


4 


00031690009 


1.62 


55600.5450 


55600.5638 










96364-01-07-00 


3.3 


55601.3754 


55601.4217 


-10.40±0.05 


0.040 


0.515 




96364-01-08-00 


3.3 


55602.3547 


55602.4000 


-10.40+0.04 


0.040 


0.517 




96364-01-09-00 


3.0 


55603.4719 


55603.5132 


-10.54+0.03 


0.029 


0.376 




96364-01-10-00 


3.2 


55604.3171 


55604.3582 


-10.63+0.06 


0.023 


0.304 




96364-01-11-00 


3.1 


55605.4297 


55605.4710 


-10.63+0.03 


0.023 


0.302 



Notes. The three vertical blocks separated by a horizontal line refer from top to bottom to April 2010, September 2010, and February 2011, 
respectively. (a) ID is the conventional name used in our broad-band spectral analysis exploiting quasi-simultaneous observations. <b) OBS IDs as 
9xxxx-0 1 -xx-xx refer to RXTE, 0003169000x to Swift, IBIS and IEM-X to the INTEGRAL monitoring. (c) 2-30 keV absorbed flux is in units of 
ergcirT 2 s -1 , not reported for instruments operating in different bands. (d) 2-30 keV luminosity in erg/s assuming a distance of 3.3 kpc. <e) different 
snapshots with the same OBS ID. 



of the source across all different states, by assuming the con- 
tinuum model for the source broad-band emission discussed in 
Sect 13.31 A log of the observations is given in Table [2] where 
we report also the source flux and luminosity. The lower en- 
ergy threshold of the RXTE/PCA is 3 keV and it does not per- 
mit a precise measurement of the absorption column density 
in the direction of RX J0440.9+4431. Thus, we have fixed it 
to the weighted average value measured from all the available 
Swift/XRT and RXTE quasi-simultaneous observations charac- 
terised by high enough S/N to significantly constrain our pa- 
rameter of interest (the observations with a conventional ID in 
Table EJ. We find Na = (2.8 ± 0.1) x 10 21 cirT 2 . In sev- 
eral fits the value of log 10 (A) in the bmc model could only be 
poorly constrained. Where possible, we estimated a lower limit 
at 90% c.l. by using the XSPEC command steppar. For the 
lower statistic spectra we fixed it to log 10 (A) = 1.5. Dur- 
ing the observations performed when the source reached the 
highest flux (observation IDs 95418-01-01-00 and 95418-01- 



02-00 , April 2010) an exponential cut-off at high energy (i.e., 
exp(-£ , /£ , poid)) was needed to properly fit the data: we have 
found ZiFoid = (30 + 6)keV and (33 + 8)keV, respectively. For 
all the other observations, we have verified that the lower limit 
on the cut-off energy is outside the energy range covered by our 
data (~50keV), consistently with the value derived using INTE- 
GRAL/ISGRI data (see Sect.l33i The results obtained from all 
the fits are reported in Fig. [4] We include in the figure also the 
EPIC-pn data from the XMM-Newton obs ervation that caught 
RXJ0 440.9+4431 in its persistent low state (lLa Palombara et all 
120121) . Note that the bmc model also provided a fully accept- 
able fit to these data Ofe- 2 /d.o.f. = 1.1/1 83) and yielded log 10 A>1.2, 
kT = (1.02 ± 0.03) keV, a = 0.81 + 0.08, r BB =(0.73 ± 0.04) km, 
and a column density A^h = (5.4 + 0.2) x 10 21 crrT 2 (uncer- 
tainties at 90% c.l.). The absorption column density measured 
with XMM-Newton/ EPIC-PN was higher than that measured 
with Swift/XRT. We verified that fitting the PCA data with this 
higher value of A^h would not affect the results discussed below. 
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Table 3. Best fit parameters for the absorbed bmc with exponential cut- 
off model for the broad band spectral fits of the 2010 outburst (uncer- 
tainties are at 90% c.L). 



N H [l0 22 crrr 2 J 
E Foid [keV] 
kT BB [keV] 
a 

A 



April 2010 a 



September 20 10 b 



lo gio 



Nbmc 
r BB d [km] 
^Gauss [keV] 

^Gauss 

CpCA 

ClSGRI 

CjEM-X2 

FluX2-30keV G 

FluXo.l-200keV 

*LM-o.f. 



0.31+0.02 

29! 4 
1.38+0.02 
0.50+0.07 
> 2.15 

030 +0 003 
u.uju 00Q1 

2.64+0.12 

6 33 +0 04 
-0.02 

(6.6+1.3) x 10- 4 
1.55+0.02 



3.78+0.01 
3.4 

1.24/571 



27 +0 03 

1 07+0.03 

i - J/ -0.04 

0.78+0.12 

0.8!°;? 

0.019±0.002 
2.13+0.16 



0.632+0.008 
0.73+0.03 
0.80±0.02 

1.272+0.009 

2.9 
0.93/559 



Notes. <a) observation A in Tabled (b) observation S in Table |2] (c) L 39 
is the X-ray luminosity in units of 10 39 erg/s, D 10 the source distance 
in units of lOkph . (d) quantity derived from Nbmc and Wbb- <c) flux 



of the RXTE/PCA observation in units of 10 erg cm 



-1 (f) 



esti- 



mated de-absorbed flux in units of 10 erg cm s 
Swift/XRT observations. 



normalized to the 



From the plots in Fig. [4] we note that the parameter log 10 (A) 
could be only poorly constrained and a shows a not very pro- 
nounced downwards trend for increasing source flux. Con- 
versely, both the effective radius and temperature of the thermal 
emission component increase with the source flux. For the effec- 
tive radius, we find (y 2 d /d.o.f.=0.6/19): 

logio = (0.182 + 0.008) + (0.39 ± 0.02) (log I0 F_ 9 ) , (3) 

where r BB is the effective radius of the black body seed pho- 
tons expressed in km and F^ is the 2-30 keV flux expressed in 

This relation is equivalent to r BB <* 
Similarly for the black body tem- 



units of 10 9 ergcm 2 s 



?0.39±0.02 



or Fv oc r. 



2.56±0.13 



perature, we find (^ 2 ed /d.o.f.=0.6/19): 

log 10 kT keW = (0.128 ± 0.004) + (0.070 ± 0.008) (log 10 F_ 



.) , (4) 



where A:7\ e v is the black body temperature expressed in keV. In 
eq. (|3]|4]i uncertainties are given at l<x c.l. To prove the varia- 
tions in the source spectral properties in the different observa- 
tions in a model-independent way, we also report in Fig. [5] the 
hardness -intensity diagram of the RXTE/PCA data (as shown 
by, e.g.. iReig & Nespolill2013l) . From this figure we note that 
during all outbursts the source X-ray emission becomes harder 
at higher fluxes and there is an evident change of the overall 
slope at PCU2 count rates >100 cts/s (the persistent emission 
level of RX J0440.9+4431 corresponds to a PCU2 count rate of 
about 1-3 cts/s) We verified that different choices of the energy 
bands would not significantly affect this result. 

An interpretation of the spectral changes in 
RX J0440.9+443 1 is provided in Sect. 03] 
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(6-10 keV)/(3-6 keV) 

Fig. 5. Hardness-intensity diagram of all the available PCA ob- 
servation of RXJ0440.9+4431. Different symbols refer to different 
outbursts: diamonds are for the outburst in April 2010 (MJD 55 292- 
55 301), stars for the outburst of September 2010 (MJD 55 447-55 450), 
and black dots for the outburst of February 201 1 (MJD 55 594-55 606). 



4. Discussion 

4. 1 . The long-term X-ray variability and the orbital period 

The Swift/BAT data allowed us to measure with a fairly good 
accuracy the orbital period of the source owing to its enhanced 
X-ray activity in 2010 and 2011. Before this period, the avail- 
able data revealed that RX J0440.9+443 1 was a persistent X-ray 
emitter at a level of about 1 mCrab in hard X-rays (15-150keV) 
and no significant orbital modulation could be detected from 
its light curve. Be/X-ray systems are rarely persistent X-ray 
sources, but many of them showed in the past only sporadic 
episodes of enhanced X-ray activity. The latter is thought to 
be originated by varia tions in the extension of the Be star equa- 
torial disk (see, e.g.. lReigll201 ll) . The existence of persistent 
Be systems supports the idea that the fainter long-term X-ray 
activity of these sources is most likely related to the accretion 
of the NS from the tenuous wind of the companion, while the 
outbursts occur as a consequence of the interaction at the system 
periastron betw een the compact o bject and the Be equatorial disk 
(see also, Rothschil d et alj|2012l , for the quiescent emission of 
A 0535+26). From the shaded regions of Fig. [T] we note that 
the maxima of the second and third flares do not correspond 
exactly to a multiple of our best candidate orbital period: one 
lies slightly before, the other after. In Be/X-ray binary systems, 
some variability of the orbital phase occurrence of outburst is ex- 
pected, but in absence of orbital ephemeris based on pulse period 
modulation, we cannot speculate further on this subject. 



4.2. The broad-band spectral model 

Previous works on persistent Be pulsars dReig & Rochel[T999l) 
showed that their spectra can usually be well described by 
a model comprising a power law component and a black 
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Fig. 4. Results of the spectral fits to the RXTE/PCA observations (filled circles) performed during the outbursts of RXJ0440. 9+4431 occurred 
in 2010 and 2011 (Table|2](. The best fit spectral parameters obtained from the XMM-Newton observation (2011 March 18) are also included 
and represented with a star (uncertainties are at 68% c.l. and lower limits at 90% c.l.).The bmc model was used for all fits, but during the two 
observations performed when the source reached the highest flux, an exponential cut-off at high energy was needed to properly fit the data (see 
Sect. l3.4t . The black body radius Tbb is obtained from eq. ([2} setting C = 1. The red solid lines in the right panels represents the best fit relations 
of Eq. l[3}|4j obtained from the PCA data only. 



body with an emitting radius o f a few hundreds meters 
(La Palombara & Mereghettil 120061 l2007t lLa Palombara et al l 
2009, 120121) . These sources are characterized by a spin period of 
the order of hundreds seconds, orbital period of hundreds of days 
and a persistent luminosity of 10 34 ~ 35 erg s _1 . The direct black 
body emission is commonly interpreted as being due to the ther- 
mal mound at the base of the neutron star accretion column (with 
a size proportional to the width of the accretion stream), while 
the power law is the signature of Comptonization of the ther- 
mal photons by the in falling plasma. Following this interpre- 
tation and using data from several X-ray observatories, we have 
verified that also the broad-band spectrum of RX J0440.9+4431 
could be reasonably well described by this phenomenological 
model, although our wide energy range requirers the introduc- 
tion of an exponential roll-over at high energies. 

A more insightful description of the parameters regulating 
the X-ray emission from the source can be obtained by using 
the bmc model. This permits to describe in a self-consistent 
way the process in which part of the thermal seed photons are 
Comptonized and up-scattered to higher energies, giving rise to 
a power-law shaped emission with part of the photons escaping 
the Comptonizzation to reach directly the observer. The model is 
a general solution of the process of Comptonization of a spher- 
ically symmetric in-falling plasma but has the major limitation 
that the high energy roll-over, due to either the thermal motion 
of electrons or the deviation from a free-fall velocity profile of 



the moving plasma, is not taken into account. This can be empir- 
ically introduce d by means of an ad ditional exponential cut-off 
(expf-E/EFniA lFarinelli et al.l2013l) . still preserving the advan- 
tage of having a self-consistent treatment of the Compton up- 
scattering process. 

Several other models describe the process of Compton up- 
scattering of seed photons, with different degrees of complica- 
tion and various hm itafions. The models comptt and compst 
(see iTitarchukl 11994 and references therein) describe a ther- 
mal Comptonization in a static plasma so, despite the possibility 
of obtaining statistically satisfactory results in terms of X-ray 
spectral fitting, they are not suitable to be applied to accreting 
sources at highly sub-Eddington rates as RXJ0440. 9+4431 (in 
these cases the accreting material is expected to be nearly in free- 
fall). 

The model comptb dFarinelli et al.ll2008l) has been proposed 
to describe the emission of low-mass X-ray binaries for both the 
cases of simply thermal or thermal plus bulk Comptonization. 
However, the correlation between the electron temperature kT e , 
the spectral index a and the bulk parameter 5 does not yield an 
independent determination of the whole parameter set. Assum- 
ing that bulk Comptonization is an important effect, the model 
implementation requires that the parameter kT e is fixed or al- 
lowed to vary in a very narrow range, inferred from other cons id- 
erations. It is also important to note that the comptt, compps, and 
comptb models do not consider the presence of a magnetic field 
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so that the classical Thomson or Klein-Nishina electron scatter- 
ing cross-sections are assumed. The Be X-ray binary nature of 
RX J 0440.9+4431, to gether with its relatively high pulsed frac- 
tion ( Usui et al.ll2012l) . suggests that the NS in this system is 
endowed with a magnetic field significantly higher than in the 
low-mass systems. This is expected to funnel efficiently the ac- 
creting material onto the neutron star's polar caps, leading to an 
electron temperature that is har dly predictable a priori. 

The model proposed by iBecker & Wolffl d2007l) is also 
not appropriate to describe the X-ray emission observed from 
RXJ0440. 9+4431, as it assumes a velocity profile for the in- 
flowing plasma in which an intense radiation field originates 
a collision-less shock above the NS surface. In this case, 
the radiation creeps through the lateral walls of the accretion 
stream giving rise to the so-called "fan beam". In the case of 
RX J0440.9+4431, the emission pattern from the NS is expected 
to assume rather the shape of a "pencil beam", as at low lumi- 
nosity the radiation is believed to be able to escape directly in 
the direction of the accretion stream and no dramatic changes 
in the pulse profile could be revealed oyer two orders of magni - 
tude in luminosity (lLa Palombara et aUl2012t [Usui et al. 2012). 
The latter assumption can be roughly demonstrated in the present 
case by comparing the magnitude of the gravitational force of the 
neutron star and the radiation pressure from the polar cap where 
matter is thermalized. We find that 



Table 4. Best fit parameters of the compmag model using a free-fall ve- 
locity profile of the accreting matter. 



1.3 x 10 14 mr 6 



(5) 



where m and r^ are the neutron mass in solar masses and the 
distance from the surface in units of 10 6 cm, respectively. On 
the other hand, the radiative force per unit mass is given by 



Prod* 1.4xlO I3 Mb b 4 /(fl), 



(6) 



where kT^B is the temperature of the thermal mound and f(B) < 
1 is a correction factor to the Thomson cross-section due to the 
presence of the magnetic field. For photons propagating along 
the vertical direction, thus almost parallel to the magnetic field, 
the Thomson cross-sec tion is strongly reduced, with f(B) <K 1 
dBecker & Wolffll2007l) . For a typical neutron star mass and us- 
ing the best fit values of kT^B reported in Tables [3] and |4] it is 

evident that & ra d <g & ar cm. 

The m odel oflBecker & Wolffl (120071) was extended recently 
(compmag; iFarinelli et al. 20121) . In this new implementation, it 
is possible to tune the velocity profile of the in-falling material. 
A broad-band spectral coverage is however necessary to obtain 
suitable constraints in the physical parameters of compmag, and 
thus the single PCA spectra could not be used for this purpose. 
More satisfactory results were obtained in the multi-instrument 
broad band spectral analysis of the outbursts in April 2010 and 
September 2010. We have chosen to use a free-fall like profile 
with fixed parameter rj = 0.5, we have also fixed the albedo at 
the NS surface to one, while the other parameters are determined 
by the spectral fit. They are the flow velocity at the NS surface 
(Bq, expressed in unities of c), the accretion column radius (ro, 
expressed in Schwarzschild radii), the seed black body temper- 
ature (kTsB), the electron temperature in the accretion stream 
(kT e ), their opticsl depth (r), and the model normalization (N). 
A summary of these results is reported in Table 2] The tempera- 
ture of the seed thermal photons is in fair agreement with the bmc 
model. The increase of the electron temperature at decreasing 
luminosity can be identified as the origin of the increased cut- 
off energy, while the normalization, N, of compmag is defined 
such that /?km = N 1 ' 2 Diq, and thus gives an estimate of the ap- 
parent blackbody emitting area. It is worth noting that compmag 





April 2010 a 


September 2010 b 


N H [l0 22 cm" 2 J 


0.36±0.01 


0.28+0.02 


kT BB [keV] 


1.50+0.01 


1.32+0.01 


kT e [keV] 


13±1 


62±1 


T 


0.12± 0.01 


0.09+ 0.01 


Bo 


> 0.56 


0.44+0.01 


>•(> 


30 +0 05 


32 +0 01 


N 


112+5 


388+1 


^Gauss [keV] 


6.33+0.03 




^Gauss 


(6.59+1.2) x 10- 4 




FluXo.i-200keV C 


3.4 


3.0 




1.21/569 


0.93/558 



Notes. <a) observation A in Table [2] (b) observation S in Table [2] (c) es- 
timate of de-absorbed flux in units of 1CT 9 ergcirT 2 s -1 , referred to 
Swift/XRT observation. 



assumes a uniform blackbody photon field distribution across 
the accreting column, so the inferred value R^, n does not actu- 
ally represent the hot spot dimension on the neutron star surface, 
but rather the total area obtained by emission at different heights 
above the neutron star surface (so that in general Ry, n > ro). This 
is confirmed here: we obtain ~ 6.5 and R^,, ~ 3.5 for the 
first and second observation, respectively (at a source distance 
of 3.3 kpc), while the column radius remains constant within the 
uncertainties at about 1 km. 

We note that the absorbing column as measured by XMM- 
Newton/PN is significantly higher than the Swift/XRT determi- 
nations, while the spectral parameters of the bmc slightly deviate 
from the trend derived from the PCA data set (Fig. This 
can be due to the different energy ranges over which spectra are 
measured which can lead to a difference in the determination 
of the spectral slope and thus of the absorbing co l umn. Fur- 
thermore, we note that according to Rom ano et al.l (120051) . the 
PN detector yields spectral parameter for the calibration source 
PSRB0540-69 which deviate from EPIC/MOS and Swift/XRT 
in the same direction as our data, i.e., higher Nh, although by 
a smaller amount. Therefore, we cannot exclude that the ap- 
parent discrepancy we have found is influenced by the calibra- 
tion systematics. We have verified that usin g a different absorp- 
tion model (tbabs) with the abundances of [Wilms et alJ ([2000) 
does not remove the discrepancy and yields an equally accept- 
able fit with compatible parameter value, except the absorbing 
column, which is systematically higher than in the model used 
in this paper. We remark that the absolute value of the absorb- 
ing column is of marginal importance for the results discussed 
here, and that the relatively low-resolution spectra analyzed here 
do not allow us to distinguish between the models. The val- 
ues we have derived a re sim ilar with the one determined by 
iTsvgankov eTail (120121) . while lUsui et~ai](l2012l) use just RXTE 
data (above 3 keV) and get absorbing column roughly one or- 
der of magnitude in excess: we have verified that excluding the 
Swift/XRT data, we would have obtained similar values, which 
we conclude to be inaccurate as they are hampered by the lack 
of low-energy coverage. 
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Fig. 6. A selection of unfolded RXTE/PCA spectra at different lumi- 
nosity states based on the bmc model. The corresponding observation 
ID is also indicated. 



Among all these different spectral models, the bmc is able to 
provide a self-consistent and robust physical description of the 
accretion process in RXJ0440. 9+4431 over the different lumi- 
nosity levels even exploiting a limited energy band. The comp- 
mag model is able to provide complementary information on 
the accretion process, which are not present in the bmc model, 
although with more structured assumptions on the geometrical 
properties of the accretion stream, in particular a significant ver- 
tical extent of the emitting plasma above the NS surface. How- 
ever, it generally requires a large energy coverage and it is of 
limited use when data are not available to provide constraints on 
the source emission above ~20keV. 



4.3. Spectral parameters 

The observations of RX J0440.9+4431 in 2010 and 2011 al- 
lowed us to study spectral changes in the X-ray emission from 
the source over two orders of magnitude in X-ray flux and mass 
accretion rates (~ 10 14-16 g s~'). To compare the results of our 
spectral analysis at different accretion rates, we plot in Fig. |6]a 
selection of our bmc best-fit unfolded spectral models obtained 
at increasing flux levels from the fits to the RXTE data. We note 
that the spectra of the lower persistent state are rather soft, they 
become harder at intermediate states, and display a more pro- 
nounced curvature at the highest luminosities. This is in agree- 
ment with the measurement of a high-energy exponential cut-off 
at ~30keV in the RXTE/PCA data during the two observations 
performed when the source flux was > 10 _9 ergcm _2 s _1 . The 
photon index of the Comptonization Green's function slightly 
decreases in comparison to the observations performed at the in- 
termediate flux level, also indicating a hardening of the source 
X-ray emission. At the intermediate state, INTEGRAL /ISGRI 
data were available to provide coverage at the higher energies, 
an exponential cut-off was measured at (53^[[j)keV. This sug- 
gests that the cut-off energy increases at lower luminosities so 
that the hard X-ray spectrum presents the hard tail expected 
wh en bulk Comptonizati on is the dominant reprocessing chan- 
nel dBecker & W olff 2005). Observations performed at the lower 



fluxes (< 10 erg cm s ) showed a steeper Comptonization 
index and a cooler thermal photon temperature, they are thus 



characterized by a significantly softer emission. These results 
are confirmed by the model-independent analysis of the source 
emission provided by the HID in Fig. [5] The observations in 
which a cut-off energy of 30keV was measured are represented 
with the two blue diamonds located in the upper right corner of 
the plot. The PCA observations performed at intermediate lev- 
els, for which quasi-simultaneous INTEGRAL /ISGRI data were 
also available, are represented with the two green stars located 
immediately below the blue diamonds. The source emission 
hardens remarkabley at count rates ~80-100cts/s, i.e., where the 
change in the cut-off energy is also revealed by the spectral fit- 
ting. The flux corresponding to the lowering of the cut-off energy 
is 1.3 t 1.8x 10~ 9 erg cirT 2 s _1 , and translates into a luminosity 
~ 2 x lO^ergs- 1 . 



According to the discussion of B ecker et al.l (1201 2|) . at low 
luminosity (Lx < 10 34 ~ 35 erg s~') matter strikes the NS surface 
after passing through a gas-mediated shock. At higher luminos- 
ity, the radiation pressure begins to affect the accretion process 
so that a radiative shock forms, and plasma deceleration takes 
place through Coulomb interactions in an vertically extended at- 
mosphere. At this stage, the efficiency of the thermal cooling 
of the high-energy electrons responsible for the Comptonized 
emission increases significantly and can produce the observed 
change in the cut-off energy. However, the luminosity level at 
which this transition happens is not well predicted on theoretical 
basis and is here tentatively identified on observational grounds 
at ~ 2 x 10 36 ergs~'. It is worth noticing that in all observa- 
tions, RXJ0440. 9+4431 remains relatively faint in the X-rays 
compared to the giant outbursts of transient Be/X-ray binaries 
(Lx ~ 10 38 ergs _1 ). As a consequence, the HID does not show 
the typical upper branch turning towards softer colors typically 
observed during the brightest outbursts of these sources, which 
might be interpreted as the onset of the fully radiation domi nated 
accretion column (Reig & Nespoli 20 131: [Becker et al]|2012l) . 

To test further our understanding of the accretion processes 
in RXJ0440. 9+4431 given above, we have estimated from the 
bmc best-fit parameters the terminal velocity of the accretion 
flow at the NS surface for all the RXTE spectra and reported 
our results in Fig. ED F or this purpose, we follow eq. (91-92) of 
iBecker & Wolff (2007) adapted to the notation used in the previ- 
ous sections of our paper. By equating the energy flux dissipated 
by the thermally emitting surface to the flux of kinetic energy, 
we obtain 



t4 ^ i 2 



(7) 



where <x SB is the Stefan-Boltzmann constant, T^b the black body 
temperature, u the stream velocity, and 



J = 



M 



(8) 



The radius of the hot spot emitting the thermal radiation, tbb, can 
be derived from the measured spectral parameters using eq. ©. 
We use also the standard conversion between mass accretion rate 
and X-ray luminosity: 



M = 



A nD 2 F x Rl 
GM t 



(9) 



where D is the source distance, Fx is approximated by F2-30keV> 
G is the gravitation constant, /?» and M» are the NS radius and 
mass. By taking into account that Nbmc — L^g/D^ (see eq. (O) 
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Fig. 7. Luminosity dependency of the velocity of the accreted matter 
at the NS surface obtained from eq.UOIassuming C = 0.25. 



and expressing the free-fall velocity at the NS surface as vg — 
^2GM t /R t , we obtain that 



v 

u ff 



N 



BMC 



4.8 x 10 7 Cf\ 



(10) 



where C is the geometrical factor of equation (|2) and Fx is ex- 
pressed in ergcm _2 s _1 . By assuming that the thermal emitting 
surface is a circular slab (C = 0.25), we derive a terminal veloc- 
ity of the accretion stream which is * 0.85 vg for all the observa- 
tions, apart for the two corresponding to the highest luminosities 
(Fig. |7]i. In these cases the terminal velocity marginally exceeds 
the free-fall limit. This might indicate that the thermal mound 
is building-up on the NS surface so that the geometrical factor 
C of the thermally emitting region increases and approaches that 
of the spherical geometry case (C=l). This indication agrees 
with the interpretation that the spectral change at the peak of 
the flaring activity of RX J0440.9+4431 is due to the formation 
of a radiatively dominated accretion column suggested by the 
decreasing cut-off energy. We note that the above calculations 
agree fairly well with the base velocity of the compmag best-fit 
parameters (J3o ~ 0.5 ~ 0.8uff/c). 

Throughout the whole luminosity range, the variation 
of the parameters of the thermal emision component in 
RXJ0440.9+4431 could be satisfactorily described by the re- 



lations: r BB <* Fg ±[lm , kT BB oc /ru.u/uiu-uos (seg Sect |J3] and 
FigHJi. In the case of accretion onto a neutron star from the stel- 
lar wind, it is expected that the radius of the hot spot is inversely 
correlated to the luminosity, as at higher accretion rate matter 
should penetrate along magnetic field lines which are more and 
more concentrated towards the NS magnetic poles dArons & Leal 
1980). Conversely, if the transfer of matter is mediated by an ac- 
cretion disk, the area of the hot spot is proportional to the ratio 

#*/#Aifven'* 4 /? ( see ' e.g- lLamb et all 19731 IWMte et al.lll983l ). 
The correlation determined from the fits to the data thus support 
the idea that RXJ0440. 9+4431 is continuously accreting from 
a surrounding disk (despite the opposite would be expected for 
the persistent activity of a Be/X-ray binary, see Sect. 14. 11 1. How- 
ever, the measured relation between tbb and the source luminos- 
ity is much steeper than expected in the case of disk accretiorQ. 
This might indicate either that matter penetrates through closed 
magnetic field lines at the border of the magnetosphere or that 



8 To infer the thermal component radius from the model parameters, 
we have assumed a constant geometry (parameter C of eq. 



the neutron star magnetic field structure is a more complicated 
than a simple dipole close to the surface. The former idea was 
suggested to justify the magnetosphe ric instabilities obse rved 
in the outbursts of EXO 2030+37 5 dKlochkov et alJl201 lh and 
A 0535 +262 dPostnov et ai]l2008l) . The latter has been invoked, 
e.g., bv lFerrigno et al.l (12009b to interpret the soft component in 
the broad-band spectrum of 4U01 15+63 and can be related also 
to accretion-induced deformations of the magnetic field and in- 
stabilities in the accretion mound which are shown to exist on the 
base of recent numerical simulations (Mukherie e et al.ll2013h . 

The way in which the accreting matter penetrates the NS 
magnetosphere and accretes onto the compact star surface is also 
producing significant changes in its spin period. In Sect. 13.21 we 
have shown that the spin period values measured from the obser- 
vations of RX J0440. 9+443 1 are not precise enough to constrain 
the source orbital properties and/or evaluate accretion torques. 
However, we note that RXJ0440. 9+4431 showe d a remark able 
spin-d own over the past 12 years from 202 s (Rei g & Rochel 
1999) to *206s. If we assume that RX J0440.9+443 1 is ac- 
creting from a disk at a non negligible rate, then the mea- 
sured spin down can be interpreted in terms of coupling be- 
tween the neutron star magnetic field lines and the external re- 
gions of the ac cretion disk located beyond the so-called coro- 
tational radius dWangj [19951 119961) . This assumption has been 
someti mes used to const rain the NS magnetic field. As pointed 
out bv lBozzo et al.l d2009l) . the uncertainties on the location and 
the extension of the boundary layer in the disk accretion mod- 
els for magnetized neutron stars are relatively large. There- 
fore, we could just derive that the magnetic field for the NS 
in RX J0440.9+4431 is of the order of fewxl0 12 G. This is in 
agreement wit h the timing signatures of the broad-band noise 
discussed in Tsyg ankov et alj d2012l) . The non-detection of 
the previously reported cyclotron scattering absorption line dis- 
cussed in Sect. 13.31 is not a counter argument to this magnetic 
field estimate, as these features are not ubiquitous in the spectra 
of accreting magnetized neutron stars. 



5. Conclusions 

In this work, we reported on the analysis of the last outbursts 
detected from the persistent Be X-ray binary RX J0440.9+4431, 
and used all archival observations of this source performed in 
2010 and 2011 to investigate its spectral properties in a broad 
X-ray energy domain (0.3-150keV). 

We obtained the most accurate measurement of the source 
orbital period so far at (150.0 + 0.2) days and performed a de- 
tailed analysis of the source spectrum covering more than two 
orders of magnitude in luminos ity. We used the bu lk mo- 
tion Comptonization model (bmc, iTitarchuk et al.llT997h . which 
provided a self-consistent picture of the accretion processes in 
RXJ0440. 9+4431. In particular, we found a coherent depen- 
dency of the size and temperature of the thermal emitting re- 
gion as a function of the luminosity. These quantities do not 
follow the relations expected for accretion onto strongly mag- 
netized neutron stars, and we argue that either accreting matter 
penetrates through closed magnetic field lines at the border of 
the magnetosphere or that the neutron star magnetic field config- 
uration deviates from a dipole close to the surface. 

Analyzing the spectral changes at the hard X-ray band 
(>10 keV), we have identified the luminosity level at which 
the accretion flow dynamics begin to be affected by radiation 
pressure (Lx ~ 2 x 10 36 ergs~'). We estimated the in-falling 
plasma velocity close to the neutron star, which is consistent 
with predictions from a physically motivated model. This pro- 
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vided an independent indication of the building-up of a ra- 
diative accretion column at the highest luminosity reached by 
RX J0440.9+4431 from the likely bending of the thermal emis- 
sion surface. RXJ0440. 9+4431 was particularly well suited 
for this study because the variations of its thermal emission 
could be unambiguously revealed during the higher and lower 
X-ray active states. This is not the case for other Be X-ray 
binaries (see, e.g., the controversy on the spectral modeling 
of EXO 2020+375 , iReig & Rochelll999t iKlochkov et al] 2007: 
IWilson et al.ll2008h . 

We compared our findings on the accretion process in 
RXJ0440. 9+4431 by exploring the most recent spectral mod- 
els accounting for thermal emission and its Comptonization in 
a strong magnetic field. The compmag model could provide a 
number of further details for the interaction between the neutron 
star hosted in this source and the in-falling plasma, but the us- 
ability of this model is limited by the statistics of the available 
data. More sensitive missions will hopefully cover the required 
broad energy range (2-50keV) (e.g., LOFT iFeroci et al.ll2012l) 
and will allow us to study in greater detail the accretion process 
at the lower end of the luminosity range of Be/X-ray binaries by 
means of complex physically motivated models, which cannot 
be significantly constrained at present. 
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